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^ ■ Abstract 

Basing on invariant properties of universal multifractals we propose a simple 



algorithm for interpolation of multifractal densities. This algorithm admits general- 



^ ' izations to a multidimensional case. Analytically obtained are multifractal charac- 



teristic of the function interpolating initial data. We establish the relation between 
the parameter existing in algorithm and the Levy index a which is the main index 



^ . for scaling function of universal multifractals 

1 Introduction 

Construction of the models for complex highly variable natural processes is closely asso- 
ciated with the study of their multifractal characteristics. Popular methods of averaging 
allows one to study smoothed characteristics only and lead to the loss of information 
about variability of the processes on scales less than the scales of averaging. Therefore, 
the problem is how to construct the function averaged over the scales e << eo in the limit 
e ^ with the help of smoothed experimental data (that is, averaged over some given 
scale eo ). In fact, it implies the reconstruction of variability at small scales. Obviously, 
this problem has no unique solution without additional requirements and principles of 
such an interpolation. Examples of fractal interpolations with different principles are 
presented, for example, in Refs. [|T|, 0. In the paper we use other principles which are 
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related to the possibility to return from interpolating function to the smoothed one by 
averaging over scales. At the same time, returning from an arbitrary interpolation step 
(that is, from the interpolating function on an arbitrarily small scale) we should get an 
initial function. The second principle is associated with the idea of similarity in a broad 
sense, and implies the conservation of memory about small-scale variability of unknown 
strict function in smoothed data. Similar principles are already widely used in various 
cascade models of turbulence (see, e.g.a-model 0, ||^, and /?-model 0, [|| ). Besides, 
we use a natural requirement that interpolated (smoothed) function possesses definite 
multifractal characteristics. This requirement is general for all methods of fractal 
and multifractal interpolation. We use the momentum scaling function K{q) as a mul- 
tifractal characteristic. It is related by the known way with the spectrum of singularity 
/(a) 01) ID) 01 and the mass power At least, we require the simplicity of the 

interpolation procedure, this requirement being hardly defined. We propose a simple and 
fast interpolation algorithm which allows one to make its generalization on a multidimen- 
sional case. The multifractal characteristics of the function interpolating initial data are 
obtained analytically. It is shown, that multifractals arising as a result of interpolation, 
belongs to universal multifractals []T2| associated with the deep properties of stochastic 
processes. The proposed interpolation method is based on the existence of the invariants 
(that is, the characteristics independent of scales) obtained in the paper for the universal 
multifractals. 



2 Formulation of the problem and the algorithm 

Let us consider the process of smoothing the function p{x). We assume without loss 
of generality that the unit interval [0,1] is the definition region of this function. Then, 
p{x) is assumed as integrable. We divide the definition region in to N equal segments 
[0,1] = U^=i -^n , In = {{n — ^)^,ne) where e = 1/A^. We defined the step-like function 
Pe{x) which has a non-zero constant value at each interval J„, 

(n+l)e 

/ p{x)dx 

P,(X) = , Xein (1) 

Where x G [0, 1] . The function Pe{x) smoothes the function p{x) at scale e . Successive 
increasing the scales of smoothing allows us pass to more and more smoothed description 
of p{x). The function pi{x) (L=l) is the last in the hierarchy of smoothed descriptions, 

L 

J p{x)dx 

Piix) = '—^ ^ ip) (2) 
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where (...) the parentheses imply spatial averaging. The function p^, a; G [0, 1] form a 
completely ordered set. It is easily proved that (...) is invariant at all stages of smoothing. 
Indeed, 



L JV-1 ('i+l)c 

/ pe{x)dx I p{x)dxe 



_ i=0 ie 



L Le 



N-l 1 L 




^ / p{x)dx = - p{x)dx = (p) 



(3) 



Therefore, the mean value of any smoothed function p^X^) coincides with that of 
the initial function. The described process of smoothing simulated the main features of 
experimental data treatment. Let us formulate the problem. We assume that the initial 
density p{x) is a multifractal with given momentum scaling function K{q). Knowing 
smoothed description Peo{x) , we try to restore Pe{x) at e << Cq with the same momentum 
function. In other words, we create a model of the functionp(a;) with conserved properties 



function of the smoothed density p^^ . According to the meaning of smoothing, it is not 
necessary for p{x) to get the fixed values which coincide with the corresponding values 
Peo{x) of at the definite x. We expect that this problem has no unique solution, because 
the integration procedure in Eq.(l) implies the loss of information about the variability 
of p{x) on scales less than e . Besides, after the reconstruction of its variability, the 
interchanges of particular segments of the function are assumed, those do not influence 
their multifractal characteristics. Therefore, let us formulate the particular algorithm 
of interpolation and discuss its properties. We assume that given smoothed step-like 
function pi (x) (p is some integer) has some value not equal zero at each interval - long. 
The absence of zero values is not principal and is assumed for the simplicity of subsequent 
analysis only. Therefore, we know Pe(a;) at e > ^, and, in particular, pi — (p) . Below we 
assume (p) without loss of generahty. We describe the algorithm as successive iterations 
of pi{x) , each of them consisting in 3 simple operations: 



1. The contraction of from the interval [0,1] onto the intervals [0,^], ,^], [^^^ ,1] 
and the construction of the function with the period ^. The periodic function repeats 

on its period the behavior of pi (a;) at the whole interval [0,1], that is it coincides with 

p 

pi{xp^) (operation P). 



2. The periodic function obtained at the 1-st step is raised to the v -th power (operation 
R). In general case the power z/ can depend on the number n of the interval = [^^, ^] 
, and on the iteration step. We restrict ourselves to the case i/=const, for simplicity. 




p 



p 
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3. The mean value is normalized (operation N), that is, the mean value at each 

interval /„ is fitted to the corresponding value of pi{x) at the same interval. This is done 

by multiplication of the values at the interval /„ on the value of pi{x) at x e /„, and 

p 

subsequent division of the function, constructed after the steps 1 and 2, on the mean value 
at this interval. 

A single iteration step of the function pi (x) consists in applying operation NRP to 

p 

it. An analogous transition is realized from the iteration step at the scale of smoothing 
efc to k+1 - st step at the scale e(^k+i), that is, 

Pe,+,(x)=NRPp,,(x) (4) 

It is seen, that during such an algorithm ek+i = e| while the number of intervals Nk+i = 
is increased according to A^^+i = N^. It implies that, if we have p intervals initially, 
then on the k-th step N/g = p^'' , that is, the variability of the function Pe^.{x) grows hyper 
exponentially. The operations P,R are based on scale similarity properties of p{x) and 
conservation of "memory" during the smoothing. Besides, R has a control parameter 
v influencing the multifractal properties of the interpolating ftinction. The operation 
N ensures conservation of the mean values and is related to the invariant of smoothing 
described above. The formulated algorithm is exclusively simple, easily realizable on 
computers and allows one to provide a complete analytical analysis of the interpolating 
functions. 

3 Analysis of multifractal properties of interpolat- 
ing functions 

Let us proceed with the study of the multifractal characteristics of the function Pe{x), 
which are constructed according to the algorithm presented above. It is convenient to 
pass from p^^ {x) to the functions defined on integers. We give a natural ordinal number 
n to the segments of the unit segment subdivision with length e^, and associate the value 
of a step-like function Pe^i^) to the function pk{n) at this segment, that is, 

pk{n) = pek{x\x e[{n-l)ek.,nek]) (5) 

Then, the index k implies the number of iteration, while n=l,2,...,iVfc. Here Nk is a 
number of the unit segment subdivision intervals with the length at the k-th iteration 
step. We formulate in an explicit form the rule of the transition (NPR) from pfc(m) to 
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Pk+i{n) {m=l,2,...,Nk; n=l,2,...,A^fe+i)- It is easily proved that 

^^''(n mod A^fc) , r n n 

Here the standard notations are used; the square brackets denote an integer part of a 
number, (j + p Nk ) mod Nk =j if, < j < A^fc while p is an integer ;n=l,..,iVfc. Mi,{k) is 
a V -th moment of Pek{^)-i that is. 



n-l n=l 

Eq. (6) allows us to calculate easily arbitrary moments of the introduced density 
Pfc(n). As an example, we prove the conservation of the mean density on an arbitrary 
iteration step. By definition, the mean density at (k+l)-st iteration step is 

^ Nk+i 



M,{k + 1) =< Pk+i >= E Pfc+lH 

Inserting the expression for pk+i{n) see Eq. (6), and taking into account that A^^+i = 
N^, after partial summation we get 

1 E Pk^in) Nk nf (u^ 1 Nk 

< >= —N.^^ZMn) = ,,(„) > 

Therefore, when using the described algorithm, the mean value is conserved at arbi- 
trary iteration step. We assume \ pi — 1 without loss of generality. Then, we determine 
how multifractal characteristics of the density Pk{n) (that is, those of Peki^)) depend on 
the free parameter v. For this purpose it is necessary to calculate more general fractional 
order moments of the constructed function pk{n). At the first step, using Eq.(6), we 
establish the relation between the q-th moment at the (k-l-l)-st iteration step with that 
at the k-th iteration step. By definition 



1 Nk+i 

M,(A; + 1) = -— E Pfc+i'W (8) 
Using Eq. (6) and A^fe+i = A^|, after partial summation we get 



M,(k + l) = ?^^M,(k) (9) 
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Here the standard notations of the moments are used, see Eqs. (7), (8). We use the 



momentum formahsm [0, [|13| in order to determine multifractal characteristics of the 
interpolating function. The main definition of this formahsm is that of the momentum 
scahng function K{q) of multifractal density pe{x) at scale resolution e: 

< p,\x) >= e-^(^) (10) 

The function K{q) is obeyed the conditions K{0) = K{1) = which are the conse- 
quences of the definition (10) and the relation \ p i = 1. It is worthwhile to note that 
K{q) is independent of scales , and related by the known way to the mass power r(g) [p^]: 

r(g) = l-q + K{q) 

The mass power is related by a Legendre transform to the singularity function /(a), see 



for example, Ref [14]. Returning to Eq.(9), we consider it as a functional equation for the 



moments of interpolating function. We look for the solution in the form 

M,{k) = e.^^") (11) 
Then, for the scaling function K(q) we get the equation 



K{q) _ aK(y)-K{qu)-K(q) 



Using efc+i = e^, it is easy to exclude the dependence on scales and get the functional 
equation for K{q)] 



2K{q) = K{qu) - qK{u) (12) 

The absence of scale dependence in Eq.(13) implies that we really get multifractal density 
p{x) as a result of interpolation. To find the solution of Eq. (13) we make a substation 
K{q) = ip{q) and get the equation for ip{q): 

2ip{q) = v{il){qu) - ^(z/)) 

We look for the solution in the class of differentiated functions. We differentiate this 
equation over q and get the following equation for the derivative ip'{q)'- 

2^\q) = u^tp'{qu) 

The solution at g > is 

ip'iq) = Koq^-'^ 
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Integrating this equation and returning to K{q) we get ultimately 

K{q) = 

Using natural boundary conditions K{0) = K{1) = 0, we exclude an arbitrary constant 



In ly 

qKq In g + ciq, if z/ = 2 



ci: 



Kiq) = { ^-'^l (13) 

where a = Kq is an arbitrary constant. It is worthwhile to note that Eq. (13) 
coincides with the known K{q) for the so-called universal multifractals |12|, |1^. Extrac- 



tion of this class has a deep mathematical origin associated with the existence of limit 
theorems for the distributions of the sums of independent stochastic quantities |]r3|. From 



the other hand, just this class of multifractals is the main pretender to the description of 
real natural processes and objects. This fact is supported by experimental investigations 
of various natural phenomena [|1^], |jl^ . 



Therefore, we have shown that interpolating function Pe{x) is universal multifractal, 
the main characteristic a = of it being the known function of the control parameter 
u. Varying u, it is possible to construct various interpolations of experimental data with 
various multifractal characteristics. It is worthwhile to note that this analysis, actually 
points to the existence of 2-parametric family of multifractal density invariants. In fact, 
it is easy to prove the following theorem: If there exist a and b such that 



^ = const (14) 



is independent on scale at g G -R, then Pe{x) is a multifractal density of universal mul- 
tifractal with the Levy index « = The proof of this theorem can easily be obtained 
if one consider Eq. (14) as a functional equation for the momentum scaling function and 
repeat the solution procedure this equation in a way analogous to that sited above. In- 
verse theorem is also valid, its proof being trivial after substitution of K{q) for universal 
multifractals. This, the proposed scheme of multifractal density interpolation is based on 
the existence of the invariants (14) (independent of scale resolution) for universal multi- 
fractal distributions. Such invariants can be useful for the study of multifractal properties 
of the turbulence and of other strongly non equilibrium systems, because these invariants 
establish relation between the moments of different orders and determine splitting the 
moments for universal multifractals. 
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In conclusion we discuss convergence of the proposed algorithm. For this purpose it 
is convenient to rewrite Eq.(9) for the scaling function. The rewritten equation describes 
the scaling function evolution during the iteration process, 

K^+i{q) = \{K^{qy)-qK^{v)) (15) 

where Kn+i{q) is the scaling function at the (n + l)-th iteration step. With the use 
of Eq.(15) we consider the scaling function at the n-th iteration step with the initial 
condition taken into account. 

i^i(?) = -^(?"-g) + A(?) (16) 

a — \ 

The function A(g) is caused by initial deviation from the scaling function of the universal 
multifractal. Then, after n iterations we get 

It follows from Eq.(17) that for i/" = 2 

i^n(g)"-^-^(g"-g) (18) 

a — 1 

if A(g) is bounded above and below (i.e. Ci < A{q) < C2 )■ Moreover, if A(g) grows 
not faster than constg^ with q growing, where f3 < a, and /3 = 1, is also included, 
then the scaling function also converges to the scaling function of universal multifractal. 
Thus, it is proved that the moments of the interpolated multifractal converge to those 
of the universal multifractal (of course, if initial deviations do not grow rather fast with 
q drawing). It implies that the proposed algorithm can be used for interpolation of 
the universal multifractals. It is noteworthy that the simplicity of this algorithm also 
allows one to the use it effectively for studying properties of the universal multifractals by 
purely analytical methods. Therefore, the proposed algorithm allows one to interpolate 
experimental data and reconstruct multifractal interpolating functions, which belong to 
universal multifractals, with arbitrary index a. It allows one to use the algorithm for 
the description and creation of the models of various natural processes and structures, 
in particular, ecologically important distributions of technogenic admixtures precipitated 
after the transport in turbulent atmosphere. 

Results 

One - parameter family of invariants of universal multifractals has been obtained. The 
existence of such invariants implies the existence of rules allowing one to express higher 
one-point moments in terms of lower ones. 
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If for a multifractal one finds one of these invariants,this signifies that such a multifractal 
is in class of the universal multifractals with a definite Levy index. A simple algorithm 
for interpolating universal multifractals based on the invariance discovered has been pro- 
posed. The convergence of the algorithm has been proved, and velocity of convergence 
has been analytically estimated. 
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